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Measured Instantaneous Viscous Boundary Layer in Turbulent Rayleigh-Benard 

Convection 
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We report measurements of the instantaneous viscous boundary layer (BL) thickness 5v{t) in 
turbulent Rayleigh-Benard convection. It is found that 5^ {t) obtained from the measured instanta- 
neous two-dimensional velocity field exhibits intermittent fluctuations. For small values, obeys 
a lognormal distribution, whereas for large values the distribution of 5v{t) exhibits an exponential 
tail. The variation of 5y{t) with time is found to be driven by the fluctuations of the large-scale 
mean flow velocity and the local horizontal velocities close to the plate can be used as an instant 
measure of this variation. It is further found that in the present parameter range of the experiment 
the mean velocity profile measured in the laboratory frame can be brought into coincidence with 
the theoretical Prandtl-Blasius laminar BL profile, if it is resampled relative to the time-dependent 
frame of Sv{t). 

PACS numbers: 44.25. -|-f, 44.20. -|-b, 47.27.-i 



An important issue in the study of fluid dynamics is 
to determine the velocity within a very thin layer in the 
neighborhood of the plates and walls, i.e. the viscous 
boundary layer (BL) [V\. The classical Prandtl-Blasius 
BL, derived more than 100 years ago for flows over a flat 
plate, remains one of the few BLs with exact theoretical 
profile based on the equations of motion. For turbulent 
boundary layers, the well-known, empirical, logarithmic 
"law of wall" can be used to describe the shape of the 
velocity proflle near the boundary. However, there are 
situations in which the BL is neither fully turbulent nor 
strictly laminar, but fluctuates with time. How to charac- 
terize such BLs quantitatively has long been a challenge. 
Here, we present a method that expresses the measured 
boundary layer quantities in a time-dependent frame that 
fluctuates with the boundary layer thickness. We use the 
turbulent Rayleigh-Benard (RB) convection as the test 
case and show that the mean velocity profile measured 
in the laboratory frame can be brought into coincidence 
with the theoretical Prandtl-Blasius laminar BL profile, 
if it is resampled relative to the time-dependent frame of 
BL thickness. More generally, the concept and the algo- 
rithm presented in this work may serve as an example 
on how to disentangle the mixed dynamics in quantities 
from local measurements. 

The Rayleigh-Benard system has become a paradigm 
for the study of general convection phenomena occurring 
in nature [31 • There already exist a number of studies 
of the viscous BL in RB convection [S- 11[, however, an 
important issue remains unresolved. Scaling wise and in 
a time averaged sense, it has been found that the mea- 
sured BL properties are the same as those of the Prandtl- 
Blasius- type laminar BL 0, @, Whereas, the time- 
averaged velocity profile is found to differ from both the 
Prandtl-Blasius-type laminar BL and the turbulent log- 



arithmic BL [10|, especially for the region around the 
thermal BL where thermal plumes are generated. We 
note that, as the viscous BL is produced and stabilized 
by the viscous shear of the large-scale mean flow, the 
fluctuations of the large-scale velocity would cause the 
viscous BL to fluctuate as well, which would in turn 
cause a fixed measurement point to be sometimes inside 
and sometimes outside the BL Jdj]. In this respect, the 
time-averaged properties are not sufficient to reveal the 
dynamic nature of the BL. As fiuctuating BLs exist in 
various flow systems, a more effective analysis method is 
thus highly desirable to unlock the intricate flow dynam- 
ics in the vicinity of plates, which is the objective of the 
present work. 

Two water-filled rectangular cells were used in the ex- 
periments, the details of which have been described else- 
where [I4I. The length, width, and height, in unit of 



cm, of the small cell are 25 x 7 x 24 and those of the 
large cell are 81 x 20 x 76. The experiment was con- 
ducted at fixed Prandtl number Pr — 4.3. The small-cell 
measurements covered the range 10^ < Ra ^ 2 x 10^° 
of the Rayleigh number Ra and the large-cell measure- 
ments were made in the regime 2 x lO^'^ < Ra ^ 7 x 10^^. 
The technique of particle image velocimetry (PIV) was 
used to measure the velocity fields above the center of the 
bottom plate and within the vertical circulation plane of 
the large-scale mean fiow. The measuring region for the 
large (small) cell experiments has an area of 17.2 x 21.5 
(11.1 X 13.8) mm^, corresponding to 63 x 79 velocity vec- 
tors, with a spatial resolution 0.27 (0.17) mm. The large 
(small) cell measurement for each Ra lasted 4.42 (3.79) 
hours in which a total of 35000 (30000) vector maps were 
acquired at a sampling rate '^2.2 Hz. The time-averaged 
statistical properties and various scaling behaviors of the 
viscous BL measured in the small cell have been reported 
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FIG. 1: (color online) (a, c) Two examples of the instanta- 
neous velocity fields measured near the center of the bottom 
plate (Ra — 1.9 x 10^^). (b, d) The horizontal velocity pro- 
files u{z, t) corresponding to those in (a) and (c), but averaged 
over —2 mm < x < 2 mm. The magnitude of the velocity in 
(a,c) is coded in both color scale and the length of the arrows 
in unit of mm/s. The solid lines in (b,d) illustrate how S^it) 
is obtained. 



previously Here we make new analysis to the data 
and present the results alongside with those from the new 
measurements in the large cell. 

Figures [Ha) and [TJc) show two examples of the in- 
stantaneous velocity fields. It is seen that the velocity 
fields near the plate show strong fluctuations and larger 
velocity magnitude yields thinner viscous BL. To obtain 
more quantitative information, we calculate this fluctu- 
ating BL thickness from the instantaneous velocity field 
using the following procedure. To reduce data scatter, 
we first coarse grain the measured velocity field by aver- 
aging it along the j:-direction over a range of 2 mm for 
the small cell and 4 mm for the large cell. This yields 
the horizontal velocity profile z, t) at time t. As the 
profile shows only weak x-dependence, we will hereafter 
present and discuss results measured at the center of the 
plate, i.e. u{z,t) = u{0,z,t) [circles in Figs. [TJb) and 
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FIG. 2: (color online) (a) A sample time trace of instanta- 
neous viscous BL thickness &v{t), normalized by its most- 
probable value (5^''. (b) PDFs of logio((5„/5r'') and (c) those 
of 5v/5^^ measured at various Ra varying from 1.25 x 10^ 
to 6.4 X 10". The solid curve in (b) marks the lognormal 
distribution for reference. 



[IJd)]. Next, to increase the statistical accuracy, each 
profile is smoothed by using the locally weighted scat- 
terplot smoothing method denoted as Us{z,t) [solid 
curves in Figs. (Ub) and[IJd)]. A linear fitting is then 
made to Us{z, t) close to the plate. Denote the maximum 
velocity of Us{z,t) as Um(i), the instantaneous viscous 
BL thickness 6y{t) is then obtained as the distance from 
the plate at which the extrapolation of the linear fitting 
of Us(z,t) crosses Um{t). The uncertainty of Sv{t), pri- 
marily due to uncertainties of Um (t) and the linear fitting 
to Us{z, t) close to the plate, is estimated to be about 3%. 
We note that the results are robust and insensitive to the 
algorithm used to smooth the profile. 

Figure [2] shows a 200-second time trace of 6v (t) , nor- 
malized by the most-probable thickness (5^^. One sees 
that the fluctuating Sy{t) exhibits intermittent features, 
i.e., the amplitude of Sv{t) can be much larger than its 
mean or most-probable value. We find that, for all mea- 
sured values of Ra, 5™^ coincides with the viscous BL 
thickness Xy (to within 4%) [ll|, defined through the 
slope of the time-averaged velocity profile at the plates 
that have been previously studied extensively. The prob- 
abihty density functions (PDF) of logio((5„/(5™P) and 
Sv/S"^P are plotted in Figs. Wli>) and [2l^c) respectively. 
For either PDF there is reasonable collapse for all mea- 
sured Ra from both cells, which suggests a universality 
of the BL dynamics in turbulent RB system. In addition, 
it is seen that the tails of the PDFs at large (> (5.™^) 
can be described by a decaying exponential distribution, 
whereas those at small {< S^^) satisfy a lognormal 
statistics. The mechanism that leads to the different 
distributions of Sy below and above its most probable 
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FIG. 3: (color online) (a) Normalized profiles of u*{z*) for all measured Ra ranging from 1.25 x 10^ to 6.4 x lO". (b) The shape 
factors Sd/Sm of u*{z*) (red symbols) and u{z) (blue symbols) vs. Ra for the small-cell (open triangles) and large-cell (solid 
circles) data. The dashed line represents the value of 2.59 for the Prandtl-Blasius BL. (c) Comparison among several velocity 
profiles: u*(z*) (red circles), u{z) (blue triangles), and the Prandtl-Blasius profile (black line). The data were measured at Ra 
= 1.9 X 10". 



values is at present unknown and would be a challenge 
for future theoretical studies. However, it is clear that 
measurement made at a fixed position in the laboratory 
frame and near the boundary layer will sample a mixed 
statistic (one pertaining to the boundary layer and the 
other pertaining to the bulk), because that position will 
be sometimes inside and sometimes outside of the fluctu- 
ating BL. To achieve a clean separation of the two types 
of dynamics in the measured velocity, we define the time- 
dependent relative vertical position z*{t) with respect to 
Sv{t), i.e.. 



and 



(1) 



and hence z* is the rescaled distance from the plate in 
unit of BL thickness. The mean velocity profile u*{z*) 
in this time-dependent frame is then calculated by aver- 
aging all values of t) that were measured at different 
discrete time t but at the same relative position z*, i.e.. 



u*{z*) = {u{z,t)\z = z*S.,{t)). 



(2) 



The results are shown in Fig. [Sfa), normalized by the 
maximum value of each profile [u*{z*)]max- One sees ex- 
cellent collapse of all profiles for the regions inside and 
around the viscous BL (z* < 2), suggesting a univer- 
sal BL profile for all Ra measured in the two cells. For 
comparison, we also calculate the time-averaged velocity 
u{z) = {u{z,t)) for each fixed position z. We note that 
u*{z*) is more universal than u(z) in the sense that pro- 
files for different values of Ra collapse better in the time- 
dependent frame than they do in the laboratory frame. 

To characterize the shape of the u*{z*) profile more 
quantitatively, we examine its shape factor 5^/ 5m, where 



u*{z*) 



[U*{z*)]r 



-}dz* 



(3) 



{1 



-}{T4^}d^* (4) 



[u*{z*)]max [U*{z*)], 



are the displacement thickness and the momentum thick- 
ness, respectively Because of the zero mean- flow in 



the central region of a closed convection cell [12|, u*{z) 
decays after reaching its maximum value and hence the 
above integrations are evaluated only over the range from 
z* = to where u*{z*) = [u* {z*)\,-aax ■ For comparison, 
we also show the shape factors of the u{z) profiles mea- 
sured in the laboratory frame, which are represented by 
blue symbols in Fig. [3jb). It is seen that for all val- 
ues of Ra, 6d/Sm of u{z) are clearly smaller than 2.59 — 
the value of a Prandtl-Blasius laminar BL profile. This 
is because the time-averaged BL quantities obtained in 
the laboratory frame are contaminated by contributions 
from the bulk that has a different dynamics. It is also 
seen that the shape factors of u*{z*) (red symbols), mea- 
sured in the time-dependent frame, are much closer to the 
Prandtl-Blasius value. This suggests that the algorithm 
using the dynamical coordinates can effectively disentan- 
gle the mixed dynamics inside and outside the fluctuating 
BLs. To see this more clearly, we directly compare in Fig. 
El^c) the velocity profiles u*{z*) (circles) and u{z) (trian- 
gles), based on the same data set (Ra = 1.9 x 10^^) but 
measured in the time-dependent and laboratory frames 
respectively jwith the theoretical Prandtl-Blasius profile 
(solid line) [l[. Here, to make a proper comparison, the 
initial slope of the Prandtl-Blasius profile is matched to 
that of the measured profiles by adjusting the path length 
xq in the similarity parameter rj = z ^/ Umax /'^xqv (i' 
is the kinematic viscosity of water.) An excellent coinci- 
dence between u*{z*) and the Prandtl-Blasius profile is 
seen for z* < 2. Whereas the values of u{z) are clearly 
lower than those of the Prandtl-Blasius profile in the re- 
gion around the viscous BL (0.7 < z/Xy <2). 
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FIG. 4: (color online) (a) Cross-correlation function g(2:, r) 
between and u{z,t) measured ai z = 0.765™''. (b) The 

normalized peak position to as a function of z/SI^^ . 
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FIG. 5: (color online) (a) Cross-correlation function g^ir) 
between Sv{t) and Um(t). (b) The normalized peak position 
To of g„(T) as a function of Ra. 



With the measured 5y{t)^ we can now study some 
of its dynamic properties and find out what is respon- 
sible for its thickness fluctuation. We calculate the 
cross-correlation function between 5v{t) and u{z,t), i.e., 
g(z,T) = {[S.v{t)- {6v{t))][u{z,t-T) - {u{z,t))]) /asau, for 
each measuring vertical position z and for each Ra. Here, 

^5 = ViMt) - {s^tW), - ^{[u{z,t)-{u(z,tw), 

and (. . .) denotes a time average. Figure |4l[a) shows two 
examples of g{z, r) measured inside the viscous BL. The 
presence of a strong negative peak for both correlation 
functions and the fact that the peaks are located at tq = 
suggest that the thinning of the viscous BL and the in- 
creasing of u{z,t) occur simultaneously. We further note 
that To = holds for z < 1.5(5™^ for all Ra investigated 
[see, e.g.. Fig. Sfb)] and hence it is a BL property. It is 
further found that Sv{t) and u{z,t) are highly correlated 
inside and around the viscous BL for all measured Ra 
[see, e.g.. Fig. IH^a)]. This suggests that the fluctuations 
of horizontal velocity measured locally close to the plate 
(e.g. at a single point) can be used as an instantaneous 
measure of the fluctuations of the viscous BL thickness 
5v{t), whereas the measurement of 6v{t) itself would re- 
quire measuring the entire velocity profile. 

From Fig. SJ^b) it is seen that for z > 1.56™^ a positive 
To is obtained, suggesting that the zero-time delay is only 
a boundary property. This can be made more clearly by 
examining the cross-correlation function gm(T) between 
Svit) and Um{t), with the instantaneous maximum ve- 
locity as a representative of the magnitude of the mean 
wind. Figure [Sf a) shows gm(T) as a function of time lag 



T. Two features are worthy of note: (i) There exists a 
strong negative correlation between 6y and Um, which 
implies that the larger Um is, the smaller Sy becomes; (ii) 
The negative peak of gmi''') is located at a positive time 
lag To, indicating that the fluctuations of Um leads the 
variation of Sy. This is because a time delay is needed 
to transfer the momentum from the mean wind to the 
interior of the BL. Figure[5jb) shows the Ra-dependence 
of the To, normalized by the typical timescale S'^p'^/v of 
momentum transfer across the BL via viscosity. It is 
seen that to/(5™p^/z^) varies around 0.14 for Ra < 10^^ 
and decreases slightly with increasing Ra for Ra > 10^^. 
Here, the decrease of to/{S^p^/i/) at high Rayleigh num- 
bers suggests an increased contribution from Reynolds 
stress to the total shear stress inside the BL llj. Direct 
study of shear stress around the viscous BL further re- 
veal that Reynolds stress becomes more important than 
viscous stress for Ra > 10^^ and hence the BL changes 
gradually from being laminar to a more turbulent one. 

Finally, we check the scaling properties of Sv{t). We 
use 6™P as the typical viscous BL thickness. The 
best power-law fits give S^^^p/H - 0.64Re"°-''^=^° °2 for 
the small cell and S^p/H - O.eORe"''-^^^^" ''^ for the 
large cell, where Re — [u* {z*)],naxH/v is the large-scale 
Reynolds number. All these results are in line with those 
obtained using the time-averaged profiles measured in the 
laboratory frame ll| and the measured 6™P/H-Re scal- 
ing exponents are in good agreement with the theoretical 
value of —0.5 for the Prandtl-Blasius BL, further confirm- 
ing the conclusion of the Prandtl-Blasius-type laminar 
BL in RB system [lll[l3|. 
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